Comment on "Phenomenological model of 
Dynamic Nonlinear Response in RelsLxor 
Ferroelectrics" 

In the recent letter [1] Glazounov and Tagantsev (GT) 
proposed the phenomenological model of dynamic re- 
sponse of relaxors. The model is based on the well 
known Landau-Khalatnikov (LK) equation for the po- 
larization dynamics dP/dt = —{l/r])dF/dP, F{P) = 
{l/2)aP'^ + ll/4)/3P'^-PE. The new element introduced 
by the authors is the averaging of the solution of LK 
equation over the broad distribution of relaxation times 
T = rj/a which, as the authors assume, resembles the low 
frequency dispersion of the dielectric response observed 
in relaxors. The purpose of the present comment is to 
clarify the physical assumptions involved in the GT ap- 
proach and show that models based on LK equation can 
not describe properly the dynamic response in relaxors 
ferroelectrics. 

Note first that LK equation completely ignores all the 
inhomogeneities of the system both temporal and spa- 
tial, and the relaxation time r is characteristic of the 
dynamics of the long wave length component of polar- 
ization with zero wave vector. The legitimate question 
to ask in this case is what is the meaning of the aver- 
aging over r proposed by GT? The following scenario 
might be true. Let us assume that that the crystal is 
divided into macroregions within which the coefficient r 
has the same magnitude. Then the average over r just 
means the average of the polarization response for dif- 
ferent macroregions. However, the above physical pic- 
ture is very far from the real situation taking place in 
relaxors, where the existence of nanoscale dynamical in- 
homogeneities (polar clusters) are well documented. The 
concept of nanoscale polar clusters provides more nat- 
ural explanation of the broad distribution of relaxation 
times in relaxors than that implied by the GT approach. 
Indeed, each polar cluster may be characterized but its 
own value of t with the correlations between the values 
of T for different clusters being negligibly small. A for- 
malism which explicitly takes into account the dynamics 
of interacting polar clusters in relaxors has been devel- 
oped in Ref[2] and applied successfully to the dielectric 
experiments in PST and PMN. 

We will illustrate below that GT approach leads to 
the values of the dielectric response which in general is 
inconsistent with the experiments in relaxors. As follows 
from GT model, the linear dielectric permittivity may be 
written in the form 

e'{LU,T)^es{T)il-q{LO,T)) (1) 

where 1 — q{uj,T) = (1/(1 +w^r^)), and (...) denotes the 
average over r. It is apparent that the quantity g(a;, T) is 
the fraction of polar regions which are effectively frozen 
on the time scale u}~^. Indeed, if for the major part of 


polar regions lot << 1, then g — > 0, and in the opposite 
limit LOT >> 1 9 — > 1. Since relaxor ferroelectrics are 
close to ferroelectric instability they are characterized by 
a strong increase of the static (field cooled) permittivity 
es(T) with the temperature decrease. At the same time 
the frequency dependent permittivity e'{LO,T) decreases 
with temperature below a well defined temperature Tm- 
The applicability of Eq.(l) in this situation assumes that 
the increase of Cs is compensated for by an increase of the 
parameter q{LO,T). However, recent NMR experiments 
[3] which obtain direct information on the magnitude of 
q{LO,T) from measurements of the inhomogeneous NMR 
line width do not support this assumption. Indeed, it 
was found[3] that in PMN, g « 0.2 at w = lO^Hz and 
T=220K which, accordingly to Eq.(l) leads to the con- 
clusion that e'(w)/es w 0.8. However, dielectric experi- 
ments [4] performed in the same sample give a much lower 
ratio e'{Lo)/es < 0.05 for the same values of the parame- 
ters. 

This example, accompanied by the given above phys- 
ical picture of the meaning of the averaging of the so- 
lution of the LK equation, clearly shows that GT phe- 
nomenological model is inconsistent with the mechanism 
of relaxor behavior since it does not take into account 
the existence of short range polar clusters. Therefore, 
one can doubt the validity for relaxors of the obtained in 
Ref [1] expression Pg ~ Pe'^{uj)e'{3Lo) for nonlinear third 
order polarization response P^ with e'(w) given by Eq.(l). 
Note, however, that GT employed the above expession 
for Pg only for evaluating the coefficient (3 using exper- 
imental values for both e'(w) and Pg. In order to test 
how sensitive is such a procedure to the particular model 
we performed the same calculations using the interacting 
polar clusters model [2]. We obtained 

P^ ~ f3e'HLo)e'{3Lo){l - q{3u:))/{l - q{Lo)). (2) 

Since in relaxors the frequency dependence q is rather 
weak, such that (?(3w) « 9('^), one can state that the 
evaluation of the coefficient /? given by GT is quite reli- 
able, although almost with the same degree of accuracy 
one could have used the expression Pg ^ Pe'^{Lo) em- 
ployed earlier in Ref [5]. 
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